Precision Measurements on a Tunable Mott Insulator of Ultracold Atoms 
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We perform precision measurements on a Mott-insulator quantum state of ultracold atoms with tunable inter- 
actions. We probe the dependence of the superfluid-to-Mott-insulator transition on the interaction strength and 
explore the limits of the standard Bose-Hubbard model description. By tuning the on-site interaction energies to 
values comparable to the interband separation, we are able to quantitatively measure number-dependent shifts 
in the excitation spectrum caused by effective multi-body interactions. 
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The observation of the superfluid-to-Mott-insulator transi- 
tion in the context of ultracold atoms lUl 121 has triggered nu- 
merous activities both in theory and in experimental physics 
|[J|. It has become clear that ultracold gaseous systems con- 
fined to optical lattice potentials are capable of serving as 
bottom-up models for condensed matter phenomena fT-'5]. In 
addition, in view of unprecedented control and read-out ca- 
pabilities i^HEl, there is justified hope that ultracold atomic 
and molecular systems will allow the implementation of quan- 
tum simulation schemes ^ (TO). While there is tremen- 
dous progress for fermionic systems confined to optical lat- 
tices fTT", 121, most experiments have so far addressed the 
case of bosonic quantum gases and in particular the quantum 
phase transition from a superfluid to an insulating Mott state 
II2I l6l 171 [TSUTSi . In this case, as long as the interaction can 
be treated as a weak perturbation, the system is described by 
the Bose-Hubbard (BH) model llT9l . One of the merits of ul- 
tracold atomic systems is the fact that all parameters of the 
BH model can be derived in a microscopic way 1 1 1. Never- 
theless, recent theoretical ll20l - E2l and experimental |23| in- 
vestigations have shown that already for comparatively weak 
interactions corrections to the standard BH model are needed. 

In this Letter we use our capability to tune interaction ener- 
gies to values comparable to the interband spacing and thereby 
leave the range of validity for the approximations of the stan- 
dard BH model description - specifically, the restriction that 
only the lowest Bloch band in the lattice is occupied, and the 
treatment of interactions via the zero-range pseudopotential 
applied in the Bom approximation. Using a Bose-Einstein 
condensate (BEC) of Cs atoms loaded into a 3D optical lattice 
potential we first investigate the superfluid-to-Mott-insulator 
transition and its dependence on the interaction strength. We 
precisely determine the on-site interaction energies includ- 
ing effective multi-body interaction shifts, demonstrating the 
breakdown of the standard approximations. Our results show 
good agreement with treatments beyond the BH model in- 
corporating both higher bands and regularization of the pseu- 
dopotential. 

The standard BH model introduces two parameters to de- 
scribe the dynamics of ultracold atoms in an optical lattice: 
the rate J/h, which describes tunneling between neighboring 
lattice sites, and the energy U, which quantifies the interac- 



tion of atoms at the same lattice site. In the presence of weak 
external harmonic confinement the Hamiltonian reads 
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where a] (a,) are the bosonic creation (annihilation) operators 
at the /-th lattice site, ni = alai is the number operator, and 
£, denotes the on-site energy shift due to an external confine- 
ment. For small U /J the ground state at zero temperature is a 
superfluid (SF), whereas for large U /J and commensurate fill- 
ing on-site interactions inhibit tunneling and the ground state 
is the Mott insulator (MI) of exponentially localized atoms. 
These limits are connected by a quantum phase transition. 
The transition point for a homogeneous system can be cal- 
culated in a mean-field approach, giving ([//y)c=:34.8 in a 
3D cubic lattice 1191 . close to the quantum Monte-Carlo re- 
suh(f//7),«29.3El. 

In the standard BH model, U and J are usually calculated 
from lowest-band Wannier functions |1]. Two-body interac- 
tions can be described via a regularized zero-range pseudopo- 
tential, as the system is dilute, and relative momenta between 
atoms are small compared with scales determined by the range 
of the interaction potential 1261 l27l . For small values of the 
s-wave scattering length a^, this pseudopotential is then re- 
placed by a 5-function in the integrals that determine U (giv- 
ing the Born approximation result). Under these approxima- 
tions, the standard BH model was successfully used to de- 
scribe a range of experiments with ultracold atoms in optical 
lattices 10 13). However, for sufficiently strong interactions the 
approximations break down. In a simple picture, two or more 
particles localized at a single lattice site with strong interac- 
tions tend to avoid each other and the on-site wave function 
increases in width to minimize the energy, resulting in cou- 
pling to higher Bloch bands. This admixture of higher bands 
results in number-dependent shifts for the on-site interaction 
energy, and the standard BH model may be modified to re- 
produce the new bound state energies II2TI by replacing U by 
a number-dependent term U{ni) 1 20l l22l424l . However, care 
must be taken, as the replacement of the pseudo-potential with 
a 5-function is in general not valid when including higher 
bands, and instead it is necessary to use the full regularized 
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FIG. 1: (color online) a) Scattering length ag for Cs atoms in the 
lowest hyperfine ground state as a function of the magnetic field B as 
calculated in Ref. 13 II 1321 . b) Example of an integrated density pro- 
file of the BEC for 50 ms TOE. The arrows indicate the EWHM d. c) 
Center peak FWHM rf as a function of lattice depth Vb for 427 ao (cir- 
cles), 320 «() (squares), and 212a() (diamonds). The solid lines are fits 
134J from which the critical lattice depth Vq is determined, d) Crit- 
ical depth Vc as a function of as- The solid line corresponds to the 
transition points for the SE-to-MI transition given by {U /J)c = 34.8 
1191 , the shaded area indicates our uncertainty for Vb- 



EIG. 2: (color online) a) Excitation spectrum in the MI phase at 
V() = 20£r and as=212a(). The solid line is a triple gaussian fit. 
We take the resonance position as the center of the gaussian peak, 
b) Frequency of the U resonance (circles) and the 2U resonance (di- 
amonds) for various values of Vq for as=212ao. The solid lines 
are the calculated frequencies corresponding to t/(I) and 2t/(I). At 
this value for aj the transition point occurs at ~ 12£'r. c) Excitation 
spectrum in the MI regime at 20£'r and as = 320ao. The solid line 
is a five-peak gaussian fit. 



zero-range potential 0281 . Note that small modifications for U 
as a function of n, are already visible for weak interactions but 
reasonably deep lattices lfT4ll23l . 

The starting point for our experiment is a BEC without 
detectable non-condensed fraction of typically 1.0 x 10^ Cs 
atoms in the energetically lowest hyperfine ground state con- 
fined by a crossed dipole trap. Atom cooling and trapping 
follow the procedures described in Ref. 1291 [30l . The cubic 
lattice is generated by three retro-reflected laser beams at a 
wavelength of A = 1064.5 nm. With the given laser power the 
maximum lattice depth Vq is SOiiR, where ER — h-^/{2mX^) is 
the atomic recoil energy with the mass m of the Cs atom. The 
strength of interactions can be widely tuned as as depends 
strongly on magnetic field B ||3T1 |32l due to multiple Fesh- 
bach resonances as illustrated in Fig.flTa). 

We first probe the transition from the SF to the MI state 
as we vary the strength of interactions, using the standard 
interference-contrast technique |2|. We gently load the BEC 
within 400 ms into the lattice, hold the atoms for 10 ms, and 
then instantly switch off both the lattice and the external trap 
to determine the momentum distribution in a 50 ms time- 
of-flight (TOE). We determine the FWHM d of the central 
peak of the resulting interference pattern 1331 as shown in 
Fig.fltb). As expected, d shows a strong dependence on Vb 
(see Fig.flTc)). We relate the onset of the MI phase to the 



abrupt kink in the data, corresponding to a critical depth Vc 
|34|. Fig.fTfd) shows Vc as a function of as- We find re- 
markable agreement with the values for the calculated MI- 
transition points |T9l for the case of a homogeneous system 
with integer density of one atom per lattice site. We note that 
in general our data on the set of transition points exhibits some 
dependence on the initial density. Here, the density is chosen 
such that for a given interaction strength the mean atom num- 
ber per lattice site in the central region of the trap is near unity 
at the transition point. 

A second signature for the MI phase is the opening up of 
a gap and hence the appearance of distinct resonances in the 
excitation spectrum |l2l[l3][35][36l. Here, we will find that the 
experiment deviates significantly from the results of the stan- 
dard BH model and the associated approximations. Fig.l2|a) 
shows a typical spectrum when the system is deeply in the MI 
phase. We plot the BEC fraction |37| after 150ms of ampli- 
tude modulation (AM) at 20 % of Vo and back-transfer into the 
initial dipole trap as a function of the AM frequency /m. The 
BEC fraction is a sensitive indicator for temperature changes 
and hence for the amount of energy deposited into the system. 
The spectrum in Fig.|2|a), taken for comparatively weak inter- 
actions and low atom density, shows three characteristic peaks 
around the energies t//2, U and 2U. The U/2 peak relates to 
a two-phonon transition, while the 2U peak corresponds to 
an excitation at the edge between the singly and doubly occu- 
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FIG. 3: (color online) a) Detailed structure of the resonance near 
U for different values of the initial atom density at aj = 427 ay (low 
density, no double occupancy (circles), intermediate density, some 
double occupancy (squares), high density, high double occupancy 
(diamonds), for details see text). Every data point is the statistical 
average of five measurements, the error bars are the standard devia- 
tion, b) Remaining atom number iij^ after AM at the same interaction 
strength as in a), normalized to wa without AM. The solid line is a 
gaussian fit. c) Remaining molecule number wm after AM at the same 
interaction strength as in a), normalized to iim without AM. The solid 
line is a double gaussian fit. 



pied shells 11351 [361 . Interestingly, when we take spectra like 
the one in Fig.|2]^a), the peaks, in particular the one at U, are 
typically not well fit by symmetric gaussian functions. The 
reason for this will become evident below. Fig.l2](b) plots the 
positions of the U and 2 U resonances as a function of Vb for 
fls =212flo and compares them to the results of the standard 
BH model. In general, the agreement is not satisfactory. For 
comparatively deep lattices (above ISZsr) we measure a sig- 
nificant downshift for both the U and the 2U resonance. Near 
the transition point (here at Vo = 12£'r) the 2U resonance is 
clearly upshifted. The latter can be understood when taking 
the spatially separated coexistence of the SF and MI phase 
into account, as for the SF part the maximal excitation prob- 
ability lies energetically above the expected value for U and 
hence also for 2U |35 36 1. In the following, we will find an 
explanation for the downshift. 

Fig.l2|c) shows an excitation spectrum as we increase the 
effect of interactions (fls=320flo). Evidently, the resonance 
corresponding to U splits into two clearly visible but not yet 
fully resolved components. The resonance at U/2 develops a 
shoulder on the high-frequency side. The splitting becomes 
more pronounced for even higher values of as- A detailed 
excitation spectrum around U is shown in Fig.[3|a) for three 
different values for the initial density (as =427 ao)- The two 



FIG. 4: (color online) a) Frequency of the upper (R2, circles) and 
lower (Ri, squares) resonance around U and of the resonance around 
2 U (diamonds) as a function of as for a lattice depth Vq = 20£r. The 
dotted lines are the calculated t/(l) and 2f/(l) values from the stan- 
dard BH model. The solid lines are the result of the more elaborate 
calculations for (/ (2) and 3(7(3) — 2(7 (2) (see text). The dashed line 
is the calculated 3 1/(3) — 2 f/(2) using the renormalized perturbation 
theory 1221 . b) Same as in a), but for Vq = 25iiR. 



components are now well separated from each other. Their 
strength depends in opposite ways on the initial density: The 
resonance at lower frequency (Ri) decreases in strength when 
the initial density is reduced and nearly disappears at low 
densities, whereas the resonance at higher frequency (R2) in- 
creases in strength. We interpret this behavior in the following 
way: Ri is caused by excitations in the doubly occupied Mott 
shell, whereas for R2 singly occupied sites are excited. We 
also detect an increased atom loss in conjunction with Ri as 
shown in Fig.lSlb). Evidently, this resonance corresponds to 
excitations of doubly into triply occupied sites, thereby lead- 
ing to three-body atom loss I.38J . In Fig.[3]^c) we plot the num- 
ber of atoms at doubly occupied sites measured directly by 
associating them to molecules via a Feshbach resonance, re- 
moving the unpaired atoms, and detecting the remaining frac- 
tion of molecules nu If39ll40l . We observe a decrease of hm at 
Ri and an increase at R2, in agreement with our interpretation 
for the origin of the two resonances. 

We map out the dependence of the resonances on the inter- 
action strength by varying the magnetic field B from the point 
at about B = 21 G, where the two resonances start to split, to 
B=40G, spanning the range from as ~ 200 ao to as ~ 900 ao 
II32II . As the loading of the lattice at high values for as leads 
to considerable heating and loss, we ramp into the MI phase 
at as =400ao, set Vq = 20£'r, and subsequently increase as to 
the desired value with a ramp speed of 5 G/ms. Fig.Wa) shows 
the frequencies for Ri and R2 and for the resonance at2U as 



a function of as- Computing the number dependent energies 
f/(2) and U{3) is non-trivial because of the anharmonicity of 
the lattice potential and the need to regularize the 5-function 
pseudopotential for the interactions. Values for U{2) were 
calculated in Ref. f2T\, and for our parameters are also well 
approximated by rescaling the exact result for two atoms in a 
harmonic trap |28| to correct for anharmonicity, using the ra- 
tio of the Bom approximation results for the lowest oscillator 
levels in our lattice and the harmonic trap Il4ni42ll . We plot 
this in Fig.Hla) as the solid black line and see that it agrees 
well with our data for R2. Note that an approach incorporating 
higher bands but not renormalizing the pseudopotential fails 
within our range of fls values. U{3) can be estimated using the 
renormalized perturbation theory of Ref. [22J. This result is 
plotted as a dashed line in Fig.Wa), and agrees well for small 
values of as. In order to reach larger values of as we would 
need to properly resum the perturbation expansion, which up 
to now is an open problem. Interestingly, we find that the 
function 3U{3)-2U{2) « U{2)/ (1 + l.34U{2)/{fico)), with 
hco the band gap, equivalent at second order in U{2)/{fico) to 
the perturbation result, agrees well with our experimental data 
for Ri, as shown by the solid red line in Fig.Hla). The same 
measurements and calculations for a lattice depth of 25 Er are 
shown in Fig.|4|b), and we see again good agreement between 
our calculations and the experimental data. For comparison 
we plot also the interaction energies calculated with the stan- 
dard BH model as dotted lines. Remarkably, this basic calcu- 
lation fits our data for U{2) even at intermediate interaction 
strengths, due to opposing corrections that arise from includ- 
ing higher bands and regularizing the pseudopotential. 

We have investigated the SF-to-MI quantum phase tran- 
sition and the MI phase over a large region of interaction 
strengths. For strong interactions, beyond single-band BH ef- 
fects appear, leading to a splitting of the excitation resonances 
in the MI phase. Our precise measurements amount to a care- 
ful calibration of system parameters, including the scattering 
length over a wide range. The splitting of the excitation res- 
onances can be used to manipulate the Mott shells indepen- 
dently, for example increasing the number of doubly occupied 
sites without loss due to excitation to triply occupied sites. 
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